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Abstract 

The action of the quantum mechanical volume operator, introduced in connection with a sym- 
metric representation of the three-body problem and recently recognized to play a fundamental role 
in discretized quantum gravity models, can be given as a second order difference equation which, 
by a complex phase change, we turn into a discrete Schrodinger-like equation. The introduction of 
discrete potential-like functions reveals the surprising crucial role here of hidden symmetries, first 
discovered by Regge for the quantum mechanical 6j symbols; insight is provided into the underlying 
geometric features. The spectrum and wavefunctions of the volume operator are discussed from the 
viewpoint of the Hamiltonian evolution of an elementary "quantum of space", and a transparent 
asymptotic picture emerges of the semiclassical and classical regimes. The definition of coordinates 
adapted to Regge symmetry is exploited for the construction of a novel set of discrete orthogonal 
polynomials, characterizing the oscillatory components of torsion-like modes. 

PACS numbers: 31.15-p; 03.65.Sq; 02.30.Gp; 04.60.Nc 
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I. INTRODUCTION 



Extension of familiar angular momentum theory to describe quantum dynamics as a 
function of discrete variables is required to cope e.g. (i)-with structure and reactivity in 
molecular, atomic and nuclear physics [H [2], (ii)-with use in quantum chemistry of elliptic 
coordinates and orbitals [31 H], (iii)-with spin network approaches to quantum gravity |5^6j. 
These approaches exploit progresses in understanding solvability of quantum systems, such 
as provided by dynamical symmetry algebras |2J E] • 

An elementary spin network picture is schematized in Fig. [T] alternatively to the tra- 
ditional sequential coupling of angular momenta, the volume operator K = Ji ■ J 2 X J3, 
first defined in [8], acts democratically on vectors J 1; J 2 and J3 plus a fourth one, J4, 
which closes a (not necessarily planar) quadrilateral vector diagram J 1 + J 2 + J3 + J 4 = 0. 
Matrix elements of K were computed in [9j to provide a Hermitian representation, whose 
features have been studied by many [XUHT2] (see also [13] for an approach based on Bohr- 
Sommerfeld quantization). Carbone et al. [13] gave a geometrically based (Ponzano-Regge 
[15], Schulten-Gordon [16]) WKB asymptotics. 

In this work, by a suitable complex change of phase, we transform the imaginary anti- 
symmetric representation into a real, time- independent Schrodinger equation which governs 
the Hamiltonian dynamics as a function of a discrete variable denoted I. The Hilbert space 
spanned by the eigenfunctions of the volume operator [17J is constructed combinatorially and 
geometrically, applying polygonal relationships to the two quadrilateral vector diagrams in 
Fig. [I] which are "conjugated" by a hidden symmetry discovered by Regge [18] . We analyze 
the consequences of this elusive and other symmetries, providing a perspective geometrical 
view, helpful for both the characterization of molecular spectra and torsion-like modes, and 
also for the extraction of the polynomial components of rotovibrational wavef unctions. 

II. DISCRETE SCHRODINGER EQUATION AND REGGE SYMMETRY 

(k) 

Eigenvalues k and eigenfunctions \t ' e of the volume operator are most simply obtained 
through the three-terms recursion relationship first introduced in [9j and analyzed in [14J, 
where analytical expressions of eigenvalues k and eigenvectors ^\ are given for Hilbert 
spaces up to dimension five. We follow the notation and units of reference [H], but find 

2 



it crucial to apply a change of phase = (— to obtain a real, finite-difference 
Schrodinger-like equation 

ai+i*^ + a^fi^A^f. (1) 

The $^ are the eigenf unctions of the volume operator expanded in the J12 = Ji + J 2 basis^] 

The matrix elements a? in (fil) are given in terms of geometric quantities, namely 

F&ji + 1/2; j 2 + l/2)F(l; j 3 + 1/2; j 4 + 1/2) 

= . , (2) 

v/(2£ + l)(2£-l) 

where F(A,B,CQ) = |p + B + C) (-A + B + C) (A - B + C) {A + B - (7)]§ is the 
Archimedes' ("Heron's") formula for the area of a triangle with side lengths A, B and C. 
Thus is proportional to the product of the areas of the two triangles sharing the side of 
length I and forming a quadrilateral of sides ji + 1, j 2 + 1, J3 + | and J4 + 1, the parameters 
entering in Eq. (J2j|^} Note that the latter physically correspond to four quantum numbers 
associated to the quantum angular momenta Ji, J2, J3 and J4. They appear to be all on the 
same footing, indicating that the volume operator can be thought of as acting democratically 
on either a composite system of four objects with vanishing total angular momentum, or a 
system of three objects with total angular momentum J 4 . 

The similarities between the discrete Schrodinger equation and the three-term recursion 
occurring for the 6j symbol motivated the authors of |14j to carry out the analysis of its 
semiclassical behavior along the lines of [T5l [T6] . Another similarity with the case of the 
6j symbol appears concerning the range of t. As noted in |17j . the basic requirement that 
the four vectors form a (not necessarily planar) quadrilateral leads to identify the range of 
i with 

D = 2 • min (j h j 2 , js, U, j'1,3'2, f s , fd + 1 (3) 

which is also the dimension of the Hilbert space where the volume operator acts. Here ji 
and j[ are conjugated by Regge symmetry (Fig. [TJ, i.e. connected by j[ = s — ji, where 

s = (ji + 32 + h + h) /2 = (j[ + 32 + fs + 3a) / 2 ( 4 ) 



1 The treatment would be analogous had we chosen the J23 = J2+J3 basis, with t = j'23. This construction 
relies on properties of the quadratic operator algebra generated by (J12) 2 , (J23) 2 and K. Once chosen 

the eigenbasis of (J12) 2 , the other two operators are also tridiagonal and have the form of Eq. ([IJ. 

2 In Loop Quantum Gravity ji labels eigenvalues of the area operator 87r L p 7-^/ (ji (ji + 1) . where L p is the 

Planck length and 7 is the Immirzi parameter (see e.g. [5j I13|). 
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FIG. 1. A quadrilateral and its Regge "conjugate" illustrating the elementary spin network represen- 
tation of the symmetric coupling scheme: each quadrilateral is dissected into two triangles sharing, 
as a common side, the diagonal I, The other sides are of length Jj = jj + 1/2 (and J[ = j\ + 1/2); 
£, which is the discrete variable in Eq. ([I]), is shown as the distance between foci 1 and 3 of the 
confocal ellipses where the vertices of the quadrilaterals lie. The two sets of four side lengths of the 
Regge conjugate quadrilaterals are obtained by reflection with respect to the common semiperime- 
ter s (Eq. Q). This relationship can be interpreted either as concerted stretchings and shortenings 
by the parameter r = (j x - j 2 + j 3 - J4)/ 2 introduced in [19], or by v = (ji - j 2 - h + J*4)/2 
occurring in the projective interpretation of Robinson |2U| . Shown is also the difference between 
the semimajor axes of the two ellipses, u = {j\ + 22 — J3 — ji)/^- Signs are decided according to 
the choice of primed and unprimed quadrilaterals. Also, u and v would exchange their roles had 



we chosen the other diagonal £ as the variable t (see note|lJ). In Eq. (11) the orthogonal nature 
of this set of transformations is exhibited explicitly by the matrix W. The passage to the Regge 
conjugate configuration (s, —u, —r. — v) is revealed as a quaternionic conjugation, motivating our 
nomenclature. 



is the semiperimeter, common to both quadrilaterals. The map between primed and un- 
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primed js is given by the symmetric 0(4) transformation in 
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(5) 



denoted i? in the following. This is a striking manifestation of the relevance of the Regge 
symmetry in the present analysis: indeed it can be checked that the volume operator is 
invariant under such symmetry and therefore its spectrum and eigenfunctions are invariant 
too. The Regge symmetry shows up also to be important to assist in determining ranges of 
i and of the js. From Eq. (|3]) one can decide to work with unprimed quantities, label as 
ji the minimum of the eight entries and ordering j 2 and j± according to ji < j 2 < 34', then 
J2 - ji = £m < t < ji + 32 + 1 = hi and j A - j 2 + ji < j 3 < U + 32 - h- The lower and 
upper limits in j'3 correspond respectively to u = and r = 0, i.e. cases when the two Regge 
conjugate quadrilaterals are coincident. 



III. HAMILTONIAN DYNAMICS 



Transparent techniques are available to study the semiclassical behavior of difference 
equations of type (JI]) (see. e.g. [2TI 122] and references therein). The Hamiltonian operator 
for the discrete Schrodinger equation ([T]) can be written, in terms of the shift operator 

H = (a e e- iip + a e+1 e lip ) with <p = -i A ( 6 ) 

representing the variable canonically conjugate to i. The two-dimensional phase space (£, (p) 
supports the corresponding classical Hamiltonian function given by 

H = 2a e+ i cos (f, (7) 

as illustrated in Fig. [2] for the two Regge conjugate quadrilaterals of Fig. [TJ now allowed to 
fold along £ with (p perceived as a torsion angle^J 

3 Similarly the dependence on i can be appreciated as a concerted bending mode, for example writing the 
product of the areas in the numerator of Eq. as (Ji J2J3J4 sin#i sin#3)/4, where 9\ and 83 are internal 
angles in 1 and 3, respectively. 
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The classical regime occurs when quantum numbers js are large and £ can be considered 
as a continuous variable. This limit for ag permits us to draw the closed curves in the (£, k) 
plane of Fig. [3] and Fig. |1J obtained when if = or tc in Eq. 0. These curves have the 
physical meaning of torsional-like potential functions 

Ut = -UJ = 2a e , (8) 

viewing the quadrilaterals in Figs. [T] and [2] as mechanical systems. Noteworthy in Fig. [3] 
is the further symmetry along the k = line ((p = vr/2), missing e.g. in the otherwise 
similar case of 6j symbol [19], where "caustic curves" are studied in a square of the (j'12, J23) 
plane (the "screen"). Here the perfect duality between j 12 and j 23 is lost, and the symmetry 
along the k = line appears in the potential functions as the continuous counterpart of the 
known fact that eigenvalues of the volume operator come in pairs, k and — k (plus possibly 
the k = eigenvalue when D is odd). Another manifestation of this symmetry links the 
eigenfunctions: it is easy to show from Eq. Q that $^ = (— 1) £ $^ , and this appears in 
Fig. [3] as the striking alternating features in the positive part of the spectrunj^J 

The phenomenology of caustics presented in [19] can be interestingly translated to this 
case, but taking into account also such additional symmetries. Remarkably Regge symmetry 
is a key to the classification of the Lissajous-type of potential functions given in Eq. 
This includes also the limits where some quantities are large, which in the 6j case lead to the 
3j's [19J. In the present case this limiting procedure can be shown to lead to the cylindrical 
or planar spin networks discussed by Neville |23j for unpolarized and polarized gravitational 
waves. 



IV. DISCRETE ORTHOGONAL POLYNOMIALS 



The preceding considerations apply to the interesting issue of extracting the polynomial 
components out of wavefunctions and to this aim the defining three-term recursion in Eq. Q 
is sufficient according to Favard theorem [24J. Actually, these polynomials can be obtained 
from the secular equation once eigenvalues are calculated. However, instead that directly 
from Eq. ([T]), we find it much more insightful (i) - to eliminate the square roots to give an 

4 The "mirror" symmetry of the Qj enlightened in |19| applies here too. In particular, allowing negative 
values of £ would permit infinite replicas of Figs. |3jand|4jon both sides of the I range. 
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FIG. 2. The two quadrilaterals of Fig. [T] looked at as a mechanical system, evolve creasing the 
pairs of triangles in which are dissected along £, according to a torsion mode corresponding to the 
same dihedral angle f + (p in both cases. Adding the edges 24 and 2' 4' two tetrahedra having the 
same volume can be visualized. In fact, their volume is proportional to H of Eq. ([T]) which is the 
product of the areas of two triangles divided by the length of the hinging edge times the sine of 
the dihedral angle. Thus classically the volume is an energy function which is a constant of motion 
along the classical trajectories which are solutions of the Hamilton equations ^ = ^ = 
Indeed, edges 24 and 2' 4' would have the same length I = j'23 had we chosen to expand the volume 
operator in the basis of J23 = J 2 + J3 (note [I]): two different confocal ellipses would describe the 
system and the vertices 2,4 would coincide with 2', 4' as the foci of the new ellipses. On the other 
hand, vertices 1 and 3 would split to give V and 3', say, lying on the new ellipses and belonging 
either to a quadrilateral or to its conjugate. 



unsymmetrical three-term recursion with polynomial coefficients; this leads at each step to a 
polynomial proportional to $^ within a normalizing factor and a phase convention; (ii) - to 
impose Regge invariance as a guideline to obtain an illuminating geometrical interpretation, 
specifically that of two triangles forming a tetrahedron, the two faces being hinged in a 
common side, and having respectively s and u, or r and v, as the other sides; (m) - to 
show that the requirement of polynomial coefficients highlights the role of the new variables 
s, u, r, v, introduced on a purely geometrical ground in Fig. [T] After some algebra, we obtain 
from Eq. (JTJ the following unsymmetrical three-term relation, which is manifestly Regge 
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FIG. 3. Two examples of spectra of the volume operator: the horizontal lines represent the eigen- 
value k, the curves are the caustics (the turning points of the semiclassical analysis), which limit 
the classically allowed region (in red U^~, in blue JJJ Eq. (pi)). As can be seen, the eigenvalues are 
symmetrically distributed with respect to the k = (which is an eigenvalue if D is odd). In green 
the stick graph of three of the eigenfunctions (unnormalized) . Left: parameters ji, j'2, js, j4=8.5, 
10.5, 13.5, 14.5 or s, u, r, v—23.5, -4.5, 1.5, 0.5. Right: all four parameters are doubled. The extrema 
of Ug~ and U7 bracket the spectrum, which can be well understood analytically and confirmed by 
extensive numerical checks. The characteristic features of Uf and U7 can be compared to those 
for the caustics for the 6j symbol |19j . 



invariant and has polynomial coefficients 

(2£ + 1) F 2 (s, u, £ - 1) pf\ + (2£ - 1) F 2 (r, v,£+l) p { e % = k (A£ 2 - l) pf\ (9) 
The relation between and the $^ of Eq. Q is given by 

(10) 



(fc) AT(h {k) A AT F(s,u,£-1) 
r f = i\£<P f and 7V^_i = — — — JSt 



F (r, v, 

a two-term relation which can be solved in closed form, see e.g. |4j. Normalization can be 
chosen by setting boundary conditions (e.g. pf^ = 1, Pg _ = 0). 



V. CONCLUDING REMARKS AND OUTLOOK 



In retrospect, we have introduced a linear transformation R which maps the two Regge 
conjugated quadrilaterals, Eq. (|5]), and another transformation W which defines the new 





FIG. 4. Potential functions U + and U~ (Eq. (|sj) ) are shown for two cases where the conjugated 
tetrahedra coincide. The cases occur when (i) - r = {i.e. j% + j'3 = j'2 + j'4) (n) - u = 
(ji + 32 = J3 + ji)- (Hi) v = (ji + j'4 = ]2 + J3). From the viewpoint of Fig. [2j case (i) 
would correspond to a "tangential" quadrilateral, while (ii) and (m) to "ex-tangential" ones. In the 
left panel, case where only one of the r, u, v variables is zero (ji, j%, j'3, j4=100.0, 110.0, 130.0, 140.0, 
v = 0), while on the right all three of them are zero. The latter is the case for an equilateral 
quadrilateral j\ = ji = js = j'4 = 120.0 (compare to [19J for the analogous cases which appear in 
the discussion of the 6j symbol). 



variables (Fig. [T| 
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(11) 



The matrix W is recognized as the famous one which provides atomic "hybrids" of tetrahedral 
symmetry by combining one s and three p hydrogenoid orbitals |25j . Also, the (s,u,r,v) 
parametrization is reminiscent of kinematic rotations [26] of the quadrilaterals whose edges 
are interpreted as distances among four equal mass bodies. 

Summarizing, the hidden Regge symmetry acts on the new variables in an interesting 
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way 



(ji,32,h,k) (fi, fz, f 3 , j'd 



(12) 



w w 



(g, u, r, v ) < ; (s,u',r'v') 
where WRW = Q = diag(l, — 1, — 1, — 1). Therefore the transformation W and the in- 
troduction of new variables permit to associate a quaternion and its conjugate to the two 
quadrilaterals twinned by Regge symmetry. Implications of this remark in the mathematical 
context of dynamical algebras will be addressed elsewhere. 

Regarding the quadrilaterals in Figs. [T] and [2] as mechanical devices, note that the u, 
v and r coordinates are interestingly analogous to parameters occurring in the Grashof 
classification of four-bar linkages [27J, the elementary moving mechanism of engines. For 
example, the conditions for identify of Regge conjugates, namely that at least one of them 
be zero, are those for the full folding of the mechanism. 

The educated guess that equivolume Regge conjugated tetrahedra are scissor congruent 
eluded a constructive proof so far [28J: our construction (Fig. [2]) works by creasing along 
a diagonal the two plane isoperimetric quadrilaterals of Fig. [TJ concertedly stretched by 
either r or u or v. The quadrilaterals are not scissor congruent: each is dissected into two 
triangles, with congruency with respect to the product, not the sum, of their areas. 



The discrete orthogonal family of polynomials of section IV is not "classical", i.e. does 
not belong to the hypergeometric families of the Askey schemes (see |29] for relevance in 
applied quantum mechanics). Suitable generalizations would be interesting to be developed: 
indeed, the situation is similar to that encountered for the Mathieu, Ince, Lame and Heun 
families occurring for the separation of variables in elliptic coordinates for the action of the 
Laplacian operator on compact manifolds, see e.g. jHED]- Limiting cases can be formulated 
accordingly, such as the cylindrical or planar spin networks [23] , or the q extensions^] 

This set of issues, besides the cases outlined here, appears to be relevant to special 
function theory, with particular reference to the development of orthogonal basis sets of 
interest in applied quantum mechanics, and specifically in atomic and molecular physics, 
and in quantum chemistry. 

The q extensions can be conveniently based on the formulation of at in Eq. (|| as a product of two 6j's 

[S] for which the q-extension is well established 
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